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I .    INTRODUCTION 

The  development  of  increasingly  sophisticated  detection 
systems  threatens  to  reduce  the  mission  effectiveness  of 
many  types  of  weapons  platforms.  Much  effort  has  been 
focused  on  methods  of  increasing  survivability  by  reducing 
detectability .  Low  observable  (LO)  platform  design  has 
therefore  become  an  important  engineering  task  over  the  last 
several  decades.  Although  the  principles  of  radar  stealth 
are  well  known,  radar  cross  section  (RCS)  prediction  and 
reduction  are  still  difficult  tasks. 

In  general,  four  basic  techniques  are  employed  in  radar 
cross  section  reduction: 

1.  shaping, 

2.  surface  material  selection, 

3.  active  cancellation,  and 

4.  passive  cancellation. 

Shaping  is  typically  the  primary  means  of  RCS 
reduction.  However,  RCS  reduction  at  one  viewing  angle  and 
frequency  is  usually  accompanied  by  an  enhancement  at  other 
viewing  angles  and  frequencies.  Surface  material  selection 
requires  the  addition  of  lossy  materials  (radar  absorbing 
material)  to  the  surface  of  the  target.  RCS  reduction  for 
this  case  is  achieved  by  the  dissipation  of  energy  inside  of 
the  material.  The  performance  of  the  platform  (e.g., 
aerodynamic    efficiency)  generally  decreases  due  to  added 


weight.  Active  cancellation  refers  to  methods  of  reproducing 
a  signal  that  is  180°  out  of  phase  from  the  target's  skin 

return  so  that  the  total  scattered  field  is  zero.  Passive 
cancellation  is  the  addition  of  a  secondary  scatterer  to 
cancel  the  reflected  field  from  the  primary  target.  Active 
and  passive  cancellation  techniques  are  usually  limited  to 
narrow  frequency  bands  and  spatial  regions. 

The  RCS  of  an  object  can  be  attributed  to  a  combination 
of  specular  and  non-specular  scattering  components.  Specular 
scatterers  refer  to  any  part  of  an  object  where  surfaces  are 
relatively  flat  in  terms  of  wavelength  and,  in  the  case  of 
monostatic  RCS,  perpendicular  to  the  incident 
electromagnetic  wave.  A  flat  and  smooth  surface  is  a  very 
strong  scatterer  at  normal  incidence,  but  the  echo  reduces 
rapidly  as  the  angle  of  incidence  moves  away  from  the 
surface  normal .  The  larger  the  surface  area  in  terms  of 
wavelength,  the  higher  the  maximum  RCS,  and  the  narrower  the 
angular  region  of  the  lobe. 

Specular  scattering  is  usually  reduced  by  shaping.  With 
proper  orientation  of  specular  surfaces,  the  scattered 
energy  can  be  deflected  away  from  the  threat  sector.  This 
leads  to  the  design  of  unconventional  platform  structures. 
In  the  case  of  aircraft,  the  design  of  a  faceted  airframe 
drastically  reduces  the  number  of  specular  directions  for 


which  large  RCS  echoes  will  be  observed.   However,  this  will 
also  create  a  multitude  of  edges. 

Edges  are  important  constituents  in  three  different 
types  of  scattering.  First,  there  is  an  edge  diffraction 
return  due  a  leading  edge  for  parallel  polarization  or  a 
trailing  edge  for  perpendicular  polarization.  Diffraction 
can  be  defined  as  any  non-specular  scattering  of  energy  due 
to  edges  and  vertices.  Next,  edges  provide  the  source  from 
which  physical  optics  end-region  return  arises.  End-region 
scattering  is  due  to  physical  optics  like  currents  at  edges 
referred  to  as  fringe  currents.  Finally,  edges  also  provide 
discontinuities  that  reflect  surface  traveling  waves, 
potentially  resulting  in  a  strong  return  in  the 
backscattering  direction.  Surface  waves  are  generated  by 
current  waves  that  flow  across  the  surface  and  are  scattered 
by  trailing  edges. 

Edge  return  becomes  important  when  high  level  specular 
scattering  is  eliminated,  which  is  the  case  for  low 
observable  platforms.  Further  RCS  reduction  may  be  possible 
by  the  geometric  control  of  edges .  A  serrated  edge  may  cause 
diffraction  and  back-scattering  to  be  dispersed  in  other 
directions,  thereby  reducing  the  RCS  in  the  threat  sector. 
Serrated  edge  treatment  may  be  applied  to  the  leading  and 
trailing  edges  of  wings.  Obvious  examples  are  the  Northrop 
B-2  and  the  Lockheed  F-117A  aircraft  shown  in  Figure  1  and 


Figure  2,   where  the  trailing-edge  zig  zags  are  clearly- 
visible  . 


Figure  1.   The  Northrop  B-2  (from  [1] ) . 


Figure  2.   The  Lockheed  F-117A  (from  [1]) 


It  would  be  desirable  if  design  guidance  based  on  the 
above  scattering  mechanisms  of  edges  could  be  provided  by  an 
analytical  approach.  However,  no  such  approach  currently 
exist.  The  objective  of  this  thesis  is  to  investigate  and 
evaluate  the  effectiveness  of  RCS  reduction  by  means  of  the 
geometric  control  of  discontinuities  with  serrated  edges. 

Chapter  II  discusses  two  major  scattering  mechanisms 
related  to  edges.  The  concepts  of  surface  traveling  waves 
and  edge  diffraction  are  explained.  Chapter  III  introduces 
three  high-frequency  techniques  commonly  used  to  compute  the 
scattering  from  metallic  edges.  They  are  the  geometrical 
theory  of  diffraction,  the  physical  theory  of  diffraction 
and  the  method  of  equivalent  currents.  Chapter  IV  contains 
a  description  of  the  software  codes  that  were  used  in  this 
thesis.  The  method  of  moments  and  its  corresponding 
implementation  in  the  computer  codes  is  also  discussed. 
Chapter  V  presents  the  data  and  the  analysis  of  the  results. 
Finally,  Chapter  VI  presents  the  conclusions  and 
recommendations  for  future  work. 


II.   NON- SPECULAR  SCATTERING 

A.     SURFACE  WAVES 

Propagation  mechanisms  that  occur  at  an  interface 
between  two  different  media  are  generally  referred  to  as 
surface  waves.  Surface  waves  are  important  in  the  absence  of 
high  level  specular  scattering,  as  in  the  case  of  low 
observable  platforms.  Surface  waves  are  non- specular ;  that 
is,  the  maximum  scattering  does  not  occur  for  the  angle  of 
observation  equal  to  the  angle  of  incidence.  Their  RCS  is 

independent  of  body  area  but  proportional  to  L/A2  ,  where  L 
is  the  length  of  the  surface  over  which  the  wave  travels  and 
A  is  the  wavelength  of  the  incident  ray  [1] . 

A  surface  wave  that  propagates  with  little  or  no 
attenuation  is  termed  as  traveling  wave.  If  the  surface  is 
terminated  with  a  discontinuity  such  as  an  edge,  the 
traveling  wave  will  be  reflected  back  toward  its  origin.  For 
an  electric  conductor,  a  traveling  wave  occurs  only  when  the 
incident  electric  field  has  component  perpendicular  to  the 
surface   and   in   the  plane   of   incidence.   The  plane   of 

incidence  is  defined  by  the  direction  of  incidence  k  and 
the  surface  normal  h  ,  as  shown  in  Figure  3 ,  where  Et  and  Hi 
are  the  incident  field  intensities.  This  polarization  is 
defined   as   transverse   magnetic   to  z  (TM  )   because   the 


magnetic  field  intensity  vector  Hi     is  transverse  to  the  z- 

axis.  If  there  is  not  electric  field  component  in  the  plane 
of  incidence,  there  will  not  be  traveling  wave  for  an 
electric  conductor. 


Figure  3.  Incident  electric  field  with  components 
perpendicular  and  parallel  to  the  surface  in  the 
plane  of  incidence  (after  [2]). 


For  a  traveling  wave,  the  illuminated  target  acts  as  a 
combination  of  antenna  and  transmission  line,  collecting 
incident  energy  and  guiding  it  along  the  surface  until  an 
edge  or  load  is  encountered.  The  edge  in  turn  generates 
backward- traveling  wave  with  significant  backscattering 
echo.  The  magnitude  of  the  reflected  echo  of  the  surface 
traveling  wave  depends  on  the  size  and  shape  of  the  edge  and 
on  the  angle  of  arrival  with  respect  to  the  edge.  In  the 
case  of  aircraft,   a  surface  traveling  wave  excites  the 


trailing  edges  of  the  wings  and  horizontal  stabilizers  when 
the  incident  waves  are  vertically  polarized.  The  leading 
edges  also  support  traveling  wave  for  horizontally  polarized 
waves  if  the  angle  of  incidence  is  within  2  0°  to  3  0°  from 
the  edge  [2]  . 

B .    DIFFRACTION 

When  a  electromagnetic  ray  impinges  on  a  finite  flat 
surface,  part  of  the  wave  is  scattered  by  the  surface  and 
part  by  the  edge  itself.  These  so-called  edge  diffracted 
rays  will  propagate  in  all  directions  perpendicular  to  the 
edge  if  the  incident  wave  is  normal  to  the  edge.  For  non- 
normal  incidence,  an  edge-diffracted  ray  can  lie  anywhere 
along  a  forward  cone  whose  half-angle  is  subtended  by  the 
edge  and  the  incident  ray.  This  is  commonly  known  as  the 
Keller  cone  after  J.  B.  Keller,  shown  in  Figure  4,  where  /3( 

and  Ps    are  the  angles  between  the  edge  and  the  incident  and 

the  scattering  rays,  respectively. 

In  contrast  to  the  ray  optics  description  of  edge 
diffraction,  it  can  also  be  viewed  as  a  disturbance  of  the 
current  near  the  leading  and  trailing  edges.  If  one  could 
determine  the  fringe  current  (i.e.,  the  change  in  current 
due  to  the  presence  of  the  edge)  ,  then  it  could  be  used  in 
the  radiation  integrals  to  compute  the  edge-scattered  field. 


This      later      view      is      a      current      based      approach      to      edge 
diffraction. 

INCIDENT  RAY 


DIFFRACTED 
RAYS 

(KELLER 
CONE) 


Figure  4.   Diffracted  ray  geometry  and  the  Keller 
cone  (after  [1] ) . 


C.    RCS  TREATMENT 

1.    Impedance  Matching 

The  reflection  and  re-radiation  due  to  surface  waves 
can  be  eliminated  from  the  threat  direction  by  matching  the 
terminating  edges.  This  is  usually  a  narrowband  approach 
based  on  the  impedance  transformation  concept  similar  to 
that  used  in  transmission  line  matching.  An  illustration  of 


10 


terminal  matching  to  remove  the  traveling  wave  is  shown  in 
Figure  5 . 


UNTREATED  EDGE 

TRAVELLING 
INCIDENT  WAVE  WAVE 

REFLECTED 

e 

GROUND  PLANE 


TREATED  EDGE 

INCIDENT  WAVE  MATCHING  SECTION 

(RESISTIVE  FILM) 

6,  " 


GROUND  PLANE 

Figure  5.   Matching  a  edge  with  a  resistive  strip 
(from  [1] )  . 


A  second  discontinuity  is  added  with  its 
characteristics  adjusted  such  that  its  scattered  field 
cancels  with  that  of  the  original  discontinuity.  Similar  to 
the  transmission  line  matching,  the  length  of  the  matching 
section  is  about  0.25A  .  The  corresponding  resistivity  value 
can  be  determined  using  the  method  of  moments  (a 
computational  technique  which  will  be  introduced  later  in 
Chapter  IV) . 

2 .    Absorber 

Another  effective  and  popular  method  of  suppressing 
surface  traveling  waves  is  to  apply  absorbers.  Although  the 
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attenuation  provided  by  these  absorbers  is  typically 
limited,  the  backward  scattering  traveling  wave  must  travel 
the  surface  twice,  and  hence  the  loss  provided  by  the 
absorbers  will  be  doubled.  The  disadvantage  of  adding 
absorbers  is  the  decrease  in  performance  of  the  platform 
(e.g.,  aerodynamic  efficiency)  due  to  additional  weight  and 
bulk. 

3 .    Geometric  Control 

In  some  cases,  it  is  possible  to  tilt  the  axis  of  an 
edge  so  that  its  high  RCS  is  out  of  the  threat  sector. 
Serration  is  a  geometric  treatment  that  can  be  applied  to 
the  leading  and  trailing  edges  of  an  aircraft's  wings. 
Examples  of  such  implementation  can  be  seen  on  the  wings  of 
the  Northrop  B-2,  as  shown  earlier  in  Figure  1.  Serrated 
edges  can  also  be  found  on  the  weapon  bay  doors  and  landing- 
gear  bay  doors  of  LO  aircraft.  The  effectiveness  of  RCS 
reduction  using  geometric  control  of  discontinuities  with 
serrated  edges  is  presented  and  analyzed  in  Chapter  V. 
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III.   APPROXIMATE  HIGH-FREQUENCY  EDGE  DIFFRACTION  TECHNIQUES 

This  chapter  presents  an  overview  of  the  three  most 
commonly  used  high-frequency  edge  diffraction  techniques. 
The  term  "high-frequency"  does  not  refer  to  the  actual 
frequency  of  the  incident  wave,  but  to  the  size  of  the 
target  when  compared  to  the  incident  wavelength. 

Geometrical  optics  (GO)  and  physical  optics  (PO)  are 
among  the  oldest  and  most  popular  high-frequency  methods 
used  for  estimating  the  RCS  of  an  object.  The  theory  of  GO 
was  based  on  the  early  studies  of  light  and  uses  simple 
formulas  involving  the  local  radii  of  curvature  at  the 
specular  point.  This  simple  approach  fails  when  the  radii  of 
curvature  become  infinite.  The  theory  of  PO  overcomes  this 
limitation  by  approximating  the  induced  surface  currents  as 
a  function  of  the  incident  field  and  integrating  them  to 
obtain  the  scattered  fields.  The  PO  approximation  holds 
provided  that  the  surface  is  relatively  large  and  the 
scattering  direction  is  not  too  far  from  the  specular 
direction.  The  PO  approximation  fails  at  wide  angles  from 
the  specular  direction  because  the  contribution  from  the 
edge  is  not  modeled  exactly. 

Over  the  years,  high-frequency  diffraction  techniques 
such  as  the  geometrical  theory  of  diffraction  (GTD) ,  the 
physical  theory  of  diffraction  (PTD)  and  the  method  of 
equivalent  currents  have  been  developed  to  calculate  the 
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Geometrical  optics  (GO)  and  physical  optics  (PO)  are 
among  the  oldest  and  most  popular  high-frequency  methods 
used  for  estimating  the  RCS  of  an  object.  The  theory  of  GO 
was  based  on  the  early  studies  of  light  and  uses  simple 
formulas  involving  the  local  radii  of  curvature  at  the 
specular  point.  This  simple  approach  fails  when  the  radii  of 
curvature  become  infinite.  The  theory  of  PO  overcomes  this 
limitation  by  approximating  the  induced  surface  currents  as 
a  function  of  the  incident  field  and  integrating  them  to 
obtain  the  scattered  fields.  The  PO  approximation  holds 
provided  that  the  surface  is  relatively  large  and  the 
scattering  direction  is  not  too  far  from  the  specular 
direction.  The  PO  approximation  fails  at  wide  angles  from 
the  specular  direction  because  the  contribution  from  the 
edge  is  not  modeled  exactly. 

Over  the  years,  high-frequency  diffraction  techniques 
such  as  the  geometrical  theory  of  diffraction  (GTD) ,  the 
physical  theory  of  diffraction  (PTD)  and  the  method  of 
equivalent  currents  have  been  developed  to  calculate  the 
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diffraction  effects  of  metallic  edges.  These  techniques  are 
applicable  if  the  edges  of  the  scattering  body  are  long 
compared  to  the  wavelength  of  the  incident  field. 

A.    GTD 

GTD  was  introduced  by  Keller  [3]  as  an  asymptotic 
method  for  the  solution  of  diffraction  problems  at  high 
frequencies.  It  was  developed  to  supplement  geometrical 
optics.  The  total  field  at  an  observation  point  can  be 
decomposed  into  GO  and  diffracted  components: 

Es=Eg0+Ed.  (1) 

The  diffracted  field  is  determined  from  a  suitable 
canonical  problem.  A  canonical  problem  is  one  that  has  the 
same  local  geometry  features  and  for  which  a  known  scattered 
field  is  available.  An  infinite  knife  is  an  appropriate 
canonical  problem  for  a  thin  plate  edge.  The  problem  of 
scattering  from  an  infinite  knife  edge  was  solved 
analytically  by  Sommerfeld  [11]  .  The  diffracted  electric 
field  can  be  expressed  as  [2], 


_rvfa 

Ed=^T—[(i'Ei)(X-Y)sx(Sxi)  +  Z0(i-Hi)(X+Y)sxi]  (2) 

sin  ft 


where  T  is  a  divergence  factor,  k-2nlX     is  the  free-space 
wave  number,  s     is  the  distance  along  the  ray  from  the  edge 
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element  to  a  far  field  observation  point,  s  is  the  unit 
vector  of  this  ray,   /?,-  is  the  angle  subtended  by  the  edge 

and  the  incident  ray  as  shown  in  Figure  4,  and  t  is  a  unit 
vector  aligned  along  the  edge.  The  diffraction  coefficients 
are  [2] , 


_     (l/n)sin(;r/n) 


Y  = 


cos(7r/«)-cos[(y/',  -y/^/n] 

(3 
(l//i)sin(7T/n) 


cos(;t/«)-cos[(i/as  +y/)/n] 

where  n  is  the  exterior  wedge  angle  normalized  with  respect 

to  K  ,   and  \f/i      and  \j/s     are  the  angles  of  the  transverse 

components  of  the  incident  and  diffracted  directions  with 
respect  to  one  of  the  surfaces  meeting  at  the  edge,  as  shown 
in  Figure  6 . 

GTD,  like  all  ray  tracing  techniques,  has  the 
limitation  of  predicting  infinite  fields  wherever  an 
infinity  of  rays  come  together  (these  points  or  lines  are 
referred  to  as  caustics) .  This  is  obviously  wrong  since  the 
far  diffracted  field  must  be  finite  in  the  real  world.  This 
problem  requires  some  form  of  "caustic  matching"  correction 
which  blends  the  wide-angle  predictions  smoothly  into  the 
known  value  of  the  caustic. 
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Region  1 
Upper  Face  Lit 


Region  3 
Lower  Face  Lit 


Figure  6.   Geometry  for  wedge  diffraction  (from  [2]). 


B.    THE  METHOD  OF  EQUIVALENT  CURRENTS 

The  axial  caustic  problem  of  GTD  can  be  eliminated  by 
the  method  of  equivalent  currents.  Note  that  any  finite 
current  distribution  produces  a  finite  result  for  the  far 
diffracted  field  when  that  distribution  is  summed  in  a 
radiation  integral.  In  addition,  the  method  of  equivalent 
currents  ensure  that  the  diffracted  fields  can  be  computed 
for  scattering  directions  out  of  the  Keller  cone. 

The  method  of  equivalent  currents  postulates  that 
fictitious  electric  and  magnetic  current   filaments   flow 
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along  the  edges .  The  magnitude  and  phase  of  the  electric  and 
the  magnetic  currents  as  proposed  by  Michaeli  are  [4] : 


=  j2(t-Ei)De      j2{t-Ht)Dem 
e      kZ0  sin2  A  it  sin2  A 


7    =       j2(t-Ht)Dm  (5 

kY0  sin  ft  sin  ft 


where  Z0=l/yo  is  the  free-space  impedance.  The  diffraction 
coefficients  De,  Dm   and  D,_  are  [2,4], 


e '         m 


—  sin—                      —  sin  — 
D  = n n + n n {g 

e  Tc-a,  w,  n-a,  yf, 

cos L-cosi-L     cos --cos— - 

n  n  n  n 


\   .   K-a,  1   .   K-a, 

sin L  .  ,        ...  x         —sin - 


m     sin  a.          tf-a,          y/.          sin  a,              n-a2  y/i 

1    cos L-cosJ-L  2         cos --cos  — 


1    .    Tt-a,                                          I    .    tt-a, 
sin L  ,  v         —sin - 


D     =     Q     rcosy/,  w  n  cos(nK-y/s)         „  n 

em     sin  ft    sin  a,          n-ax          y/           sina2              n-a2  y/ 

r»  '     cos L-cos  —  2         cos --cos  — 


n  n 


with 


(2  =  sin  ft  cot  ft  -  sin  ft  cot  ft  ( 9 
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[sin2ft-sin2ftcos2y/g1/2 

sinaj=- ^ £ Zi±—  (10) 

sin  ft 


and 


11/2 


[sin  ft  -  sin  ft  cos  (n/r  -  wt  )] 

sina2=- ^ ^s         Ys  (ll) 

sin  ft 


With  the  method  of  equivalent  currents  proposed  by 
Michaeli,  edge-diffracted  fields  remain  finite  in  caustic 
directions  and  the  scattering  direction  is  no  longer 
confined  to  a  generator  of  the  Keller  cone. 

C.    PTD 

PTD  is  a  current-based  technique  to  obtain  the  edge 
scattered  field  of  metallic  bodies.   It  was  developed  by 
Ufimtsev  based  on  a  refinement  of  the  physical  optics  (PO) 
approximation  [10]. 

Ufimtsev  postulates  that  the  total  induced  surface 
current  Js  ,     is  the  sum  of  the  PO  contribution  J p0  ,     and  an 

edge  contribution  J f  : 

Js=JP0+Jf  (12) 

Ufimtsev  refers  to  JP0     as  the  "uniform"  current  and  Jf     as 

the  "non-uniform"  current.    If  the  non-uniform  current  is 
known,   then  the  diffracted  field  due  to  the  edge  can  be 


calculated  by  integrating  this  current  over  the  surface. 
However,  the  exact  non-uniform  current  can  only  be 
determined  by  solving  a  boundary-value  problem  which  is  very 
complicated. 

In  order  to  obtain  the  edge-diffracted  field,  Ufimtsev 
subtracted  the  PO  current  from  the  total  induced  current, 
which  can  be  determined  from  a  canonical  problem.  The 
approximate  canonical  problem  for  edge  scattering  is  an 
infinite  knife  edge.  Using  analytical  methods,  Js  is 
determined  for  the  knife  edge,  then  Equation  (12)  used  to 
find  the  residue  Jf    which,  finally,  is  used  as  an  estimate 

for  the  fringe  current  on  a  finite  edge. 

Although  PTD  is  simple  in  principle,  it  can  be  rather 
difficult  to  implement  in  practice.  Complicated  two- 
dimensional  integration  is  often  necessary  as  the  fringe 
currents  can  extend  quite  a  distance  from  the  edge.  In 
addition,  a  comprehensive  theory  for  fringe  current 
interaction  has  yet  to  be  developed. 

As  a  final  comment,  it  can  be  shown  that  all  of  these 
high  frequency  diffraction  techniques  are  equivalent  [5]. 
That  is,  in  the  limit  of  increasing  frequency,  GTD  and  PTD 
give  identical  results.  An  advantage  of  these  methods  is 
that  they  provide  some  physical  insight  into  the  scattering 
from  edges  and  can  be  used  to  predict  the  direction  of 
maximum  RCS .  However,  they  are  limited  to  long  edges   (in 
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terms  of  wavelength) .  For  the  serration  geometries  of 
interest  in  this  research,  which  are  in  the  order  of  IX    or 

less,  these  methods  are  not  suitable  for  computing  low-level 
edge  diffraction.  A  rigorous  approach,  such  as  the  method  of 
moments,  should  be  used.  The  method  of  moments  technique  is 
discussed  in  the  next  chapter. 
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IV.   COMPUTER  CODES 

A.    PATCH 

PATCH  is  a  method  of  moments  (MM)  triangular  subdomain 
code  that  was  developed  by  Sandia  Labs  [6]  .  It  is  a  FORTRAN 
computer  code  that  computes  electromagnetic  scattering  and 
radiation  based  on  the  MM  solution  of  the  E-field  integral 
equation  (EFIE) .  It  models  a  continuous  surface  by  discrete 
triangular  patches  in  order  to  handle  arbitrary  shaped 
bodies  easily  and  accurately.  This  chapter  presents  a  brief 
summary  of  the  method  of  moments  and  its  implementation  in 
PATCH . 

1.    Method  of  Moments 

The  method  of  moments  reduces  the  EFIE  to  a  set  of 
linear  equations  that  can  be  solved  using  standard  matrix 
methods.  The  number  of  unknowns,  and  hence  the  size  of  the 
matrix  equation  that  must  be  solved,  depends  on  the  number 
of  triangular  patches  that  are  used  to  represent  the 
scattering  body.  In  general,  the  more  triangules,  the 
shorter  the  edge  lengths  and  the  more  accurate  the  results. 
A  useful  rule  of  thumb  is  to  limit  the  edge  lengths  to  0.1A  . 

A  surface  or  patch  based  method  of  moments  solution  is 
rigorous;  that  is,  as  the  triangles  become  smaller,  the 
method  of  moments  solution  converges  to  the  correct  value. 
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Unlike  a  wire  simulation  (wire  grid  model) ,  the  area  between 
edges  which  form  the  triangle  facets  are  solid  material,  not 
air  gaps.  Therefore,  current  truly  flows  on  the  surface  of 
the  object,  not  just  along  wires  approximating  the  surface. 
The  unknown  quantity  in  the  EFIE  is  the  surface  current 

distribution  Js.      To  apply  the  method  of  moments,  Js      is 

represented  by  a  series  of  unknown  coefficients 


^=XV„  d3) 


n=\ 


where  { J n  }  are  the  basis  functions.  Typical  basis  functions 

include  pulses,   8  functions,   triangles  and  sinusoids.  To 

determine  the  unknown  coefficients   { In  } ,   a  matrix  of  N 

linear  equations  is  formed.  By  exploiting  the  properties  of 

the  basis  functions,  the  { In }  are  determined  and  thus  Js    is 

also  determined.  Once  the  current  coefficients  have  been 
determined  using  the  method  of  moments  procedure,  the 
radiation  patterns  and  the  scattered  fields  can  be  computed. 
The  application  of  the  method  of  moments  depends  on  the 
capability  and  capacity  of  the  computer.  Because  of  the  edge 
length  restriction  for  convergence,  bodies  with  large  number 
of  triangles  generate  matrices  that  are  too  large  for  the 
computer  or  require  processing  time  that  are  too  long  for 
practical  application.  In  general,  the  number  of  edges  and 
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unknowns  (i.e.,  number  of  simultaneous  equations  that  must 
be  solved)  is  almost  the  same.  Small  differences  occur  due 
to  the  fact  that  there  are  more  than  two  current 
coefficients  associated  with  some  edges.  On  the  other  hand, 
if  the  edge  of  a  triangle  is  also  the  edge  of  the  body,  then 
there  is  no  current  flowing  across  it  and  therefore  no 
expansion  coefficient  is  required.  PATCH  uses  overlapping 

rooftop  functions  for  Jn  .   The  current  at  a  point  on  a 

triangle  is  the  vector  sum  of  the  currents  crossing  the 
edges  weighted  by  their  distances  from  the  edge. 

Several  output  files  are  generated  by  PATCH.  They 
include  the  RCS  as  a  function  of  the  spherical  angles  6  and 
0  which  are  shown  in  Figure  7 . 


P(x,y,z) 


Figure  7.  Spherical  coordinate  systems  (from  [1]) 
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2.  PATCH  Input 

The  input  file  for  PATCH  is  an  ASCII  file  that  contains 
all  of  the  geometry  information  and  calculation  parameters. 
The  file  can  be  generated  using  the  preprocessing  code  BUILD 
that  is  distributed  along  with  PATCH.  BUILD  is  capable  of 
generating  basic  geometrical  shapes  and  combining  them  to 
yield  more  complex  shapes.  In  addition  to  geometry 
information,  BUILD  also  prompts  the  user  for  calculation 
information  such  as  frequency,  observation  angles,  and 
excitation  conditions. 

The  user-input  data  is  appended  to  the  geometry  file  and 
written  in  a  format  that  is  recognizable  to  PATCH.  Upon 
execution,  the  Naval  Postgraduate  School  (NPS)  version  of 
PATCH  (which  has  been  modified  from  the  original)  looks  for 
a  file  named  "inpatch"  in  the  same  directory.  It  performs 
the  required  calculations  and  generates  an  ASCII  output  file 
named  "outpatch"  which  contains  all  of  the  output  data  in 
readable  form.  An  output  file  of  current  coefficients  is 
also  generated. 

3 .  PATCH  Code  Capabilities 

Apart  from  building  the  shape,  PATCH  provides  some 
options  to  the  user,  one  of  which  is  the  symmetry  option. 
This  option  enables  the  user  to  reduce  the  number  of 
computations  and  thus  the  time  required,  by  using  image 
planes  placed  on  the  x=0,  y=0  or  z  =  0  planes.  These  image 
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planes  can  be  either  perfect  electric  conductors  (PEC)  or 
perfect  magnetic  conductors  (PMC) .  Other  quantities  such  as 
voltage  sources  or  distributed  impedance  follow  the  image 
properties  and  are  taken  into  account  automatically.  The 
code  capabilities  are  summarized  in  Table  1. 

PATCH  also  has  the  option  of  entering  the  number  of 
impedance  loads,  faces  with  surface  impedance  loads  and 
computing  the  field  at  near  observation  points.  All  these 
options  make  the  PATCH  code  a  powerful  and  useful  tool  in 
the  calculation  of  RCS  for  a  variety  of  structures. 


Table  1.  Summary  of  PATCH  Capabilities 


Arbitrary  Shape 

Open/Closed  objects 

Modeled  by  triangular  "patches" 

Variable  patch  density 

Front  end  for  graphical  composition 
Arbitrary  edge  multiplicity 
Non-orientable  surfaces 
Symmetry  planes  may  be  included 
Multiple  bodies 

Surface 

Basis  functions  yield  surface  currents 

Type:  Wilton  Rau 

Free  of  line  and  point  charges 
Equivalent  Thevenin  circuits 
Lumped  and  surface  impedance  loading 


Excitation 

Voltage  sources 
Plane  waves 
Both 

Calculated  Quantities 

Surface  currents 
Far  field  patterns 
Radar  cross  section 
Field  calculations  at 
general  observation 
points 

Frequency  Domain 

Pattern  loops 
Frequency  loops 


B.    ACAD 

The  geometry  preprocessor  BUILD  is  only  capable  of 
simple  shapes  such  as  plates,  cylinders,  cones,  and  spheres. 
For  computational  EM  codes   to  be  practical  engineering 
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tools,   an   efficient   surface  meshing   operation  must   be 
available . 

The  Advanced  Computer  Aided  Design  (ACAD)  [7,8] 
application  is  capable  of  performing  the  desired  auto- 
meshing.  Databases  from  other  CAD  programs  can  be  imported 
into  ACAD  and  then  meshed  and  output  in  a  special  "facet" 
format.  The  facet  file  is  ASCII  and  contains  the  node  and 
facet  information  required  by  PATCH,  although  it  is  not  in 
the  proper  format.  An  ACAD-to-PATCH  translator  was  written 
to  reformat  the  facet  file  into  one  that  is  recognized  by 
PATCH.  The  relationships  between  codes  are  summarized  in 
Figure  8 . 

ACAD  provides  users  with  the  ability  to  create  and 
modify  geometry  in  two  or  three  dimensions.  Users  can  choose 
to  model  geometry  with  wireframes,  surfaces,  or  solids.  Data 
entry  to  ACAD  is  accomplished  through  one  of  many  input 
modes  available  to  the  designer.  ACAD  models  can  be  viewed 
orthographically  or  perspective.  Users  can  specify  view 
orientation  and  choose  to  display  geometry  in  multiple 
window  configurations.  Window  operations  such  as  panning, 
zooming,  and  auto  extents  are  accomplished  at  any  time 
providing  instream  capability.  The  ACAD  user  can  also 
control  the  display  of  surfaces  or  solids  with  options  as 
wireframe,  hidden  line  removal,  flat,  or  Gouraud  shading. 

At  the  heart  of  the  ACAD  system  is  the  associative 
database.  In  an  associative  database,   geometry  is  linked 
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together  in  a  relational  structure  that  remembers  parent/ 
child  dependencies.  This  type  of  database  enables  rapid 
modifications  of  geometry,  since  modifying  one  geometric 
element  automatically  adjusts  its  dependencies  based  on  a 
set  of  predefined  rules.  For  instance,  changing  an  edge  of  a 
wing  will  automatically  regenerate  any  surface (s)  built  with 
the  edge  spline.  In  turn,  all  geometry  that  is  associated 
with  the  modified  wing  (i.e.,  plane/curve  and  surface 
intersections,  fillets)  will  automatically  regenerate. 
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program 
name 

input  file 

output 

name 

file  name 

*  denotes  an  arbitrary  name  subject  to 
each  code's  character  restriction. 
Extensions  must  be  included  as  indicated. 


atp.x      - 


in.facet       *.facet 


out.  patch 


bldmat.x 


C 


outpatch 


MATLAB 


».     knit.x 


ACAD    - 


'.facet 


i 


.iges 
Mgs 


AutoCad 


mbuild.x 


i 


inpatch 


mpatch.x 


*.m  pattern 
files 


JL 


plot  patterns 


,m  geometry 
files 


plotaeometry 
1         pltpatcn.m 


in. patch 


out.facet 


pta.x 


— 
i  i 


in.  patch 


CODES 

AutoCad,  ACAD:  computer  aided 
drawing  (CAD) 

atp.x:  ACAD  to  inpatch  translator 
pta.x:  inpatch  to  ACAD  translator 
bldmat.x:  build  matlab  geometry 
files 

mbuild.x:  build  simple  geometries; 
specify  excitation  parameters  for 
mpatch.x  calculations. 

mpatch.x:  MM  current  and 
scattering  calculations 

knit.x:  remove  duplicate  edges; 
check  for  geometry  file  errors 


Figure  8.  Flow  chart  illustrating  the  relationships 
between  the  various  computer  codes  (from  [6]). 


28 


V.    RESULTS  AND  ANALYSIS 

It  is  hypothesized  that  further  reduction  of  the  RCS 
beyond  that  obtained  by  orienting  the  surface  can  be 
achieved  by  means  of  the  geometric  control  of 
discontinuities  using  serrated  edges.  Although  the  use  of 
serrated  edges  for  RCS  reduction  can  be  clearly  seen  on  the 
Northrop  B-2  aircraft  shown  in  Figure  1,  and  was  mentioned 
in  several  papers  and  references  [1,2],  not  much  data  on  the 
reduction  magnitude,  the  associated  geometry,  or  the  design 
methodology  are  available  in  the  open  literature.  Parameters 
of  interest  include  the  number  of  basic  serration  cells 
(triangles)  required  per  wavelength,  and  the  aspect  ratio  of 
the  triangles  that  form  the  zig  zags . 

In  order  to  analyze  and  verify  the  effectiveness  of 
serration  in  RCS  reduction,  an  infinitely  thin  metallic 
plate  in  the  x-y  plane  is  considered.  The  RCS  of  such  a 
plate  with  serrated  edges  is  computed  and  compared  against 
the  RCS  of  a  plate  of  the  same  size  without  serrated  edges. 

To  avoid  early  detection  by  an  enemy's  radar,  the  RCS 
contributions  from  a  LO  aircraft  at  low  grazing  angles  are 
particularly  important,  and  have  to  be  as  low  as  possible. 
Translating  this  consideration  to  a  thin  metallic  plate  in 
the  x-y      plane,   analysis   is   focused  at  angles   in  the 

vicinity  of  0  =  ±70°.  The  infinitely  thin  assumption  is 
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valid  if  the  wing  of  the  aircraft,  which  is  represented  by 
the  plate,  is  thin  compared  to  the  wavelength.  Supersonic 
aircraft  tend  to  have  thinner  wings  than  subsonic  aircraft. 

A.     SINGLE-EDGE  SERRATED  PLATE 

Analysis  was  first  performed  with  serration  applied  to 
only  one  side  of  the  infinitely  thin  metallic  plate.  The 
generic  geometry  of  such  a  plate  is  given  in  Figure  9,  where 
a  is  the  dimension  of  the  plate  along  the  x-axis,  b  is  the 
dimension  of  the  plate  along  the  y-axis,  c  is  the  dimension 
of  the  serrated  edge  along  the  y-axis,  and  M  is  the  number 
of  serration  cells  within  a  .  All  of  the  dimensions  are 
specified  in  terms  of  A  . 


0- 


X 


Figure  9.  Geometry  of  an  infinitely  thin  metallic 
plate  in  the  x  —  y  plane,  with  isosceles 
triangular  serrated  edges. 
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Based  on  the  geometry  given  in  Figure  9,  a  RCS 
reduction  of  about  10  dB  can  be  attained  for  vertical 
polarization  at  angles  around  0  =  -70°  to  -90°.  The  results 

are  shown  in  Figures  10  and  11.  The  serrations  applied  to 
the  trailing  edge  successfully  reduced  the  backscattering 
due  to  surface  traveling  waves.  However,  no  significant 
reduction  is  observed  in  the  case  of  horizontal 
polarization. 
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Figure  10.    RCS  and  reduction  level  with  a  =2 A, 
b=2l,     c  =0.25  A  and  M  =2  . 


31 


Vertical  polarization,  theta-theta 
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Figure  11.    RCS  and  reduction  level  with  a  =  X, 
b=2l  ,     c=0.5A  and  M  =1 . 


After  an  inspection  of  the  results  for  the  geometry  of 
Figure  9,  variations  in  the  serration  shape  were  examined  to 
obtain  improvement  in  the  RCS  for  horizontal  polarization. 
This  may  be  possible  if  the  number  of  edges  can  be  reduced 
without  reducing  the  number  of  serrations  M.  This  is  based 
on  the  assumption  that  the  RCS  in  the  transverse  plane  due 
to  diffractions  from  the  edges  may  be  reduced  if  the  number 
of  edges  is  reduced.   This  can  be  accomplished  by  using 
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to  diffractions  from  the  edges  may  be  reduced  if  the  number 
of  edges  is  reduced.  This  can  be  accomplished  by  using 
right-angle  triangular  edges  instead  of  isosceles  triangular 
edges,  as  shown  in  the  generic  geometry  given  in  Figure  12. 


Y 

n 


4^\ 

x 


Figure  12.  Geometry  of  an  infinitely  thin 
metallic  plate  in  the  x-y  plane,  with  right- 
angle  triangular  serrated  edges . 


Based  on  the  geometry  shown  in  Figure  12,  and  with 
reference  to  the  results  given  in  Figures  13  and  14, 
significant  RCS  reduction  in  the  case  of  horizontal 
polarization  was  achieved  with  the  serration  applied  at  the 
leading  edge.  The  amount  of  reduction  from  0  =  70°  to  90°  is 

kept  well  above  2  0  dB.  At  certain  angles,  the  reduction  can 
be  as  high  as  32  dB.  In  the  case  of  vertical  polarization, 
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reduction  level   is   greater   than  10  dB  around  0  =  -7  0°  to 
-90°,  as  shown  in  Figures  15  and  16. 


Horizontal  polarization,  phi-phi 
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Figure  13.   RCS  and  reduction  level,  where  a  =4  A, 
b  =  X,     c=0.5774A  and  M  =4  . 
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Horizontal  polarization,  phi-phi 
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Figure    14.       RCS    and    reduction    level,    where    a  =5 A, 
b  =  X,     c=0.5774A    and    M  =5  . 
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Vertical  polarization,  theta-theta 
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Figure  15.   RCS  and  reduction  level,  where  a  =4  X, 
b  =  X,     c=0.5774A  and  M  =4  . 
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Vertical  polarization,  theta-theta 
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Figure  16.   RCS  and  reduction  level,  where  a  =5  A, 
b=l,     c=0.5774A  and  M  =5  . 


In  the  case  of  the  vertical  polarization,  reduction  is 
achieved  by  applying  the  serration  at  the  trailing  edge.  The 
serrated  edges  effectively  reduced  the  RCS  contributions  due 
to  surface  traveling  waves.  As  for  the  horizontal 
polarization,  reduction  is  achieved  by  applying  the 
serration  at  the  leading  edge.  The  serrated  edges  reduced 
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the   RCS   contributions   due   to   edge   diffraction   in   the 
transverse  plane. 

In  order  to  investigate  an  overall  reduction  for  both 
vertical  and  horizontal  polarizations,  the  design  of  a 
double-edge  serrated  plate  is  considered  in  the  next  phase 
of  analysis. 

B.    DOUBLE-EDGE  SERRATED  PLATE 

With  reference  to  the  results  obtained  based  on  the 
geometries  shown  in  Figures  9  and  12,  a  double-edge  serrated 
plate  was  created  to  investigate  the  RCS  reduction  for  both 
vertical  and  horizontal  polarizations.  To  reduce  RCS 
reduction  in  the  case  of  horizontal  polarization,  right- 
angle  triangular  edges  are  applied  at  the  leading  edge.  As 
for  the  case  of  vertical  polarization,  isosceles  triangular 
edges  are  applied  at  the  trailing  edge.  The  generic  geometry 
of  such  a  double-edge  serrated  plate  is  shown  in  Figure  17 . 

RCS  reduction  for  both  polarizations  was  achieved  with 
the  geometry  given  in  Figure  17.  The  results  are  shown  in 
Figures  18  and  19.  However,  note  that  the  reduction  level 
provided  by  the  double-edge  serrated  plate,  for  either  the 
vertical  or  horizontal  polarization,  is  less  than  that 
provided  by  the  single-edge  serrated  plate.  Reduction  in 
both  polarizations  is  possible  with  the  double-edge  serrated 
plate,  whereas  the  larger  reduction  in  magnitude  for  a 
single  polarization  occurs  with  the  single-edge  serrations. 
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Hence  the  application  of  serrated  edges  is  a  compromise,  and 
it  depends  on  the  specific  mission  profile  of  the  LO 
aircraft  and  the  threat  radar's  polarization. 
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Figure  17.  Geometry  of  an  infinitely  thin 
metallic  plate  in  the  x— y  plane,  with  double- 
edge  serrations;  right-angle  triangular  serrated 
edges  applied  at  the  leading  edge  and  isosceles 
triangular  edges  applied  at  the  trailing  edge. 
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Horizontal  polarization,  phi-phi 
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Figure    18.         RCS    and    reduction    level    where     a =4  A, 
b  =  X,     c=0.5774A,     J=0.5774A,     M  =4    and    N=0.5A. 
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Vertical  polarization,  theta-theta 
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Figure  19.    RCS  and  reduction  level  where  a  =4 A, 
b  =  X,     c=0.5774A,  J=0.5774A,  M  =4  and  W=0.5A. 
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VI .   CONCLUSIONS 

The  four  basic  techniques  typically  employed  in  radar 
cross  section  reduction,  namely  shaping,  surface  material 
selection,  active  and  passive  cancellation,  have  been 
discussed  in  Chapter  I.  The  significance  of  edge  scattering 
to  the  case  of  LO  platforms  has  also  been  discussed.  The 
edge  serration  method  can  be  considered  as  a  combination  of 
shaping  and  passive  cancellation.  The  concept  of  two  major 
scattering  mechanisms  related  to  edges,  namely  surface 
traveling  waves  and  edge  diffraction,  have  been  explained. 
The  three  most  commonly  used  high-frequency  techniques  for 
predicting  scattering  from  metallic  edges  (GTD,  PTD  and  the 
method  of  equivalent  currents)  have  been  presented.  Although 
the  approximate  high-frequency  methods  were  not  used  to 
compute  the  RCS  of  the  plates  in  this  research,  they  provide 
a  method  for  predicting  the  location  and  level  of  edge- 
scattered  fields.  Futhermore,  the  ray  model  provides  some 
physical  insight  into  the  behavior  of  diffracted  fields. 
These  analytical  techniques  are  generally  used  for  edges 
that  are  5  to  10  A  or  greater  in  length,  although  they  can 
still  give  useful  results  for  edges  as  short  as  1  or  2  A  . 

The  objective  of  this  thesis  is  to  investigate  and 
evaluate  the  effectiveness  of  RCS  reduction  by  means  of  the 
geometric  control  of  discontinuities  through  serrated  edges. 
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The  results  and  analysis  have  been  presented  in  Chapter  V. 
It  has  been  known  for  sometime  that  the  use  of  relatively 
large  serrated  edges  (i.e.,  about  10  A)  is  an  effective 
means  of  RCS  reduction.  Such  application  can  be  seen  on  the 
zig  zag  trailing  edge  of  the  Northrop  B-2  and  the  Lockheed 
F-117A  aircraft.  Of  particular  interest  in  this  research  is 
the  effectiveness  of  relatively  short  edge  serrations  in  the 
reduction  of  the  principal  plane  RCS.  Because  of  the  short 
edge  lengths,  the  method  of  moments  was  used  to  compute  the 
surface  current  and  scattered  field.  The  results  obtained 
show  that  significant  RCS  reduction  for  low  grazing  angles 
is  achievable  with  smaller  serrated  edges  (i.e.,  about  lA 
per  serration  cell) . 

The  potential  of  RCS  reduction  by  means  of  the 
geometric  control  of  discontinuities  using  serrated  edges 
has  been  clearly  demonstrated.  In  practice,  the  serrations 
would  probably  be  applied  to  the  surface  of  the  wing  as 
illustrated  in  Figure  20.  Most  modern  aircraft  are 
constructed  from  composite  materials,  which  are  non- 
conducting. A  layer  of  conducting  foil  can  be  laid  over  the 
composite  structure.  The  conductor  would  also  provide 
electromagnetic  shielding  for  the  internal  aircraft  parts. 

Further  research  should  be  conducted  to  determine  the 
bandwidth  of  the  RCS  reduction.  By  performing  several 
calculations  at   frequencies  above  and  below  the  design 
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frequency,  an  effective  bandwidth  of  the  treatment  can  be 
obtained.  Finally,  the  spatial  distributions  of  the 
diffracted  fields  are  also  interesting.  This  research  only 
examined  the  diffracted  field  in  the  principal  planes.  The 
reduction  of  RCS  in  one  direction  is  generally  accompanied 
by  an  increase  in  other  directions  unless  cancellation  or 
absorption  has  occurred.  A  contour  plot  of  RCS  would  provide 
information  on  the  redistribution  of  the  scattered  field. 


Figure  20.  Serrated  edges  applied  to  the  surface  of 
an  aircraft's  wing. 
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